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Abstract. - We describe spatio-temporal correlations and heterogeneities in a kinetically 
constrained glassy model, the Kob- Andersen model. The kinetic constraints of the model alone 
induce the existence of dynamic correlation lengths, that increase as the density p increases, 
in a way compatible with a double-exponential law. We characterize in detail the trapping 
time correlation length, the cooperativity length, and the distribution of persistent clusters of 
particles. This last quantity is related to the typical size of blocked clusters that slow down 
the dynamics for a given density. 



Glassy dynamics in kinetically constrained systems has been a subject of wide interest 
during the last years [1]. It indeed includes a large class of models of particles constrained by 
dynamical rules where the relaxation is slowed down when the density of particles increases 
or the temperature decreases, sharing similarities (though with a number of differences) with 
spin-glass models or glass models [2] when the temperature is decreased. In the former class, 
slowing down of the dynamics can be understood as a consequence of steric effects, that make 
certain moves of particles impossible because of an effective high local density, putting strong 
constraints on the number of allowed paths between configurations [3]; whereas for the latter 
class, frustration and/or disorder in the interactions produces a slow time relaxation. In real 
glasses both effects are likely to be present. While most studies on the glass state have focused 
for a long time on the latter class of models, more attention is paid to the fact that the former 
mechanism is also of importance, especially since the works of Harrowell [4] on spatio-temporal 
heterogeneities. These heterogeneities are also of great interest for experimentalists, since they 
have been identified in several glassy systems [5] . 

In this paper we will focus on the numerical study of the Kob- Andersen (KA) model [6] . 
In this 3-dimensional model of lattice gas on a cubic lattice, a particle can hop to one of 
its empty nearest neighbor sites only if it has strictly less than m occupied nearest-neighbor 
sites, and if after the hop it will also have less than m occupied nearest-neighbor sites (which 
ensures detailed balance). In the original study of the KA model, the value of m is chosen 
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to be TO = 4, which is a compromise between free diffusion of particles with excluded volume 
(to = 6) and a truly non-ergodic model (for m < 3 small cubes of occupied sites are blocked 
for ever). The relaxation is hence a diffusive process with additional steric constraints. The 
energy landscape is trivial in the sense that there is no interaction between particles, hence 
all configurations are thermodynamically equally probable. However the original numerical 
study of the KA model showed features similar to those seen in glassy systems such as the 
arrest of the diffusion constant at some density pc smaller than 1 [pc — 0.881 in [6]), and 
slow stretched-exponential relaxation of the intermediate scattering function; though, it was 
already pointed out by the authors that the Mode-Coupling predictions for glasses are actually 
not valid for this model. Moreover, the dynamics is time translationally invariant (there is no 
aging), because all configurations are equivalent. One way to induce aging is to do density 
crunches (i.e quenches to higher densities) in a system where a reservoir of particles is added, 
as studied in [7], or to add explicit two-body interactions between the particles. 

Recently the question of finite-size effects in the KA model, already addressed in [6] , came 
back, related to the still unsolved issue of a true dynamical transition in the thermodynamic 
limit. It is now established that in the thermodynamic limit and in finite dimension, no 
dynamical arrest occurs at a density smaller that one [9-11]; whereas in a mean- field ver- 
sion of the model on a Bethe lattice, such a dynamical transition at pc < 1 occurs. More 
precisely, the critical density in three dimensions depends on the size L of the system like 
1 — Pc{L) = C/(ln(lnL)), where C is a constant, which makes the convergence to 1 very 
slow and unreachable by numerical simulations. In other words, at fixed density p, in order 
to eliminate finite size effects, one has to study system sizes much larger than the spacing 
between mobile particles S(p) (x exp(exp(C/(l — p))). 

Finally the KA model has been shown to share an important feature with the dynamics of 
mean- field spin-glass models and supercooled liquids [5,12-15]: the dynamical susceptibility 
x(t) shows a pronounced maximum at an intermediate time i*, this maximum is interpreted 
as the time where heterogeneities are the strongest in the system, and increases as the density 
increases (analog to a decrease of temperature in other models), hinting at the existence of 
a growing correlation length. Because of its purely constrained nature (no energy barriers, 
no dynamical transition) the origin of these spatio-temporal heterogeneities in the KA model 
comes merely from steric effects [9, 16]. 

In this study we give a description of the spatio-temporal correlations and heterogeneities 
that characterize the dynamics of the KA model. First, we characterize the distribution of 
trapping times of particles (for recent work on different systems see for example [17,18]) and its 
spatial correlations. Then we describe the behavior of the four-point spatio-temporal correla- 
tion functions in real space, giving a determination of a cooperative length-scale not provided 
a priori by x(t), in an approach complementary to [19]. Finally we study the dynamics of 
the KA model in terms of the distribution of persistent clusters (this is a new approach to 
the problem), and show how spatial heterogeneity and cooperativity can be identified in this 
way. We show that the different lengths identified in this study behave in the way predicted 
by Toninelli, Biroli and Fisher [9-11] as a function of the density. 

Trapping times. - We have computed the distribution of trapping times (or persistence 
times) P{t). Let be the time after which the particle starting from site i first moves 
to one of its nearest-neighbor site: we calculate the probability for to have the value r, 
and we average over all particles and over a large number of initial configurations of the 
system. We have computed P{t) for different densities (ranging from 0.5 to 0.8, see figurePl. 
Our fits are consistent with a stretched exponential behavior for the integrated probability 
distribution, dt P{t) = exp { — (r/rtrap)''}, where (3 decreases as the density increases 
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Fig. 1 - P{t) versus r in for different values of p. The smooth curves are for the stretched exponential 
best fits for the integrated probability. In the inset Cp(r) for p = 0.77. 



from 0.80 {p — 0.5) to 0.45 {p = 0.8). Such stretched exponential behavior also characterizes 
P(t) in KCM, in Lcnnard- Jones systems and in experiments (see [8] and references therein). 
In particular, in the East model, a temperature dependent stretching exponent for P(t) is 
observed [20]. 

The site dependent trapping times defined above do not need to be homogeneous in space. 
One can ask whether the KA dynamics induces a spatial correlation length among them. 
We have quantified this aspect by computing spatial correlations: Cp{r) = J2ai('''sTs+p) — 
(ts) {'''a+r)) : whcrc (• ■ •) stands for an average over initial configurations. lc{p) is a dynamical 
coherence length such that for distances larger than Idp) the trapping times are uncorrelated. 
An exponential fit, where Cp(r) ~ Q-^f-M ^ done for r going from 1 to 5 works well. Although 
in the density range 0.7-0.8 lc{p) is small and of the order of one or two lattice spacings, the 
growth of Icip) for increasing p (from 0.7 to 0.77) can be fitted by a double-exponential law 
l,{p) ~ 0.17exp(exp(0.16/(l - p))). 

The finding of this stretched exponential behavior leads to consider the relaxation of cor- 
relation functions. Ref. [21] has analyzed stretched exponential relaxations in glassy systems, 
looking at it as a consequence of a renewal process for a single effective particle, where the 
distribution of renewal times is itself a stretched exponential, leading to stretched-exponential 
relaxation with the same exponent for the time correlations. It is not obvious whether the 
cooperative dynamics of the KA model could be mapped onto such a simple renewal process 
with a distribution of times equal to the distribution of persistence times P{t). To test this, 
we define the overlap function as q{t) = j^j^^rp^^i{ni{t)ni{0) — p^), where ni{t) = or 1 
is the occupation number of site i at time t. We have computed {q{t)) as an average over 
initial configurations. Even if in the asymptotic time regime {q{t)) should decay as a power 
law [22], we do not see any sign of this behavior (^): it is in fact very well fitted by a stretched 
exponential, with a stretching exponent that depends on p, {q{t)) ^ exp{-~{t/treiaxY^)- 7 
decreases from 0.8 to 0.35 for p increasing from 0.5 to 0.8. The exact values of 7 are not 
identical to the ones of (3, but they are very similar: in the limits of the statistical and sys- 
tematic error, that we estimate to be smaller than 0.05, we find the values compatible. A 
possible real discrepancy of the two exponents would be a signal of the fact that the growth 
of spatial correlations with density (although lc{p) is always small in our simulations) plays 
an important role, and that the KA model can not be reduced to a unique effective renewal 
process. Instead, one has to take into account the spatial correlations of the dynamics. 



(^)We have checked that even in the unconstrained lattice gas for high density values we are not able to 
observe a simple power law decay. 
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Fig. 2 - (?4(r, t) for p = 0.85 and L = 24 as a function of time. From top to bottom, r = 1, 2, 3, 4, 5. 

Spatio-temporal correlation functions. - The dynamic susceptibility x(i) = N {{q'^{t)) — 
has been one of the main evidences presented in [15] stressing the relation of the KA model to 
liquids and spin glasses: as a function of time it shows a maximum at t* , where the sensitivity 
of the system is maximal. We have repeated with larger precision, for a range of densities, the 
measurements of [15], finding compatible results. Both the value of the maximum of and 
t* as a function of p are very well fitted by the double-exponential scaling form of [10] (see 
also Figure EJ). At t = t* (which is of the same order of magnitude as the relaxation time) 
the statistical heterogeneity of the system is maximal. We show now that it is also related to 
spatial cooperativity and heterogeneity. 

We define and measure a space-dependent susceptibihty 34 (r, t), which generahzes x(i) [14]: 

g4{r,t) ^ {Np^il - pfY' J2 {{n^{t)nMnJ{t)n,{0)) - {n,{t)n,{0)){n,{t)n,{Q))) . (1) 

\ri—rj \ —r 

Here we only take the spatial separation on a single coordinate axis, so that r is at the same 
time an Euclidean and a Manhattan distance. We show 34 (r, t) as a function of time in figure 
121 for all values of r it shows a maximum. The position of the maximum is of the same order 
of magnitude as t* but shifts to larger times as r increases. The value of the maximum is 
decreasing with r. 

The function 54 (r, t) gives information on how the decorrelation between times and t at 
site i is related to the decorrelation of site j between the same times. The decrease of (;J^"^(r) 
(in figure enables us to determine a cooperativity length above which the decorrelating 
events are not correlated (or not "cooperative"). We fit to the form g^''^{r) = exp~^^^^P\ 
The fit works well enough: for p in the range 0.7 — 0.86 we find values of a(p) and ^(p) ranging 
from 0.22 to 0.64, and from 0.70 to 2.75 respectively. Again the cooperativity length ^(p) is 
very well fitted by a double exponential form: £^{p) ~ 0.34 exp(exp 0.16/(1 — p)). We have 
determined a time scale and a cooperativity scale ^{p) that diverge as p ^ 1. Note that the 
behaviour of g™"^ (r) is more reminiscent of a strong glass of the type of the 3d Fredricksen- 
Andersen model, than of a fragile glass (see [24]), although the KA model is also fragile, but 
with a conserved order parameter. 

Connected clusters of frozen particles. - In exploring length and time scales of KA we 
have also focused on the analysis of the connected clusters of blocked particles: we reconstruct 
the clusters of nearest-neighbor particles that have not moved at t from their initial position. 
We have analyzed them by computing their number and the probability distribution of clusters 
of size n, P{n, t). We show in figure^a the number of such blocked clusters, that turns out to 
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Fig. 3 - gT'^^i'i") as a function of the distance r. 
0.86, 0.85, 0.83, 0.8, 0.77, 0.75, 0.73, 0.7. 



From top to bottom, 



be maximum at a time truster ■ We plot with different lines data at different density values. 
This time scale, close to Tf^ap and t^eiaxi has again a behavior compatible with [10]. 

At small times the probability P{n^ t) is mainly composed of two distinct populations: 
large blocked clusters that dominate the rate of relaxation and small blocked clusters. As 
time passes some particles move, and the large initial blocked clusters decrease in size, giving 
rise to smaller clusters. At an intermediate timescale (which is very close to t duster) the two 
populations of blocked clusters merge. We show in figure^lb P{n, tduster) for different density 
values. Here the behavior looks as a power law where fi ~ 1.3) and one can observe 

a collapse of the curves for all densities; the origin of this unique exponent is unclear to us. 
At very long times the system relaxes completely: the persistent clusters only consist of a few 
particles, that will eventually fly. For t tduster the initial structure of clusters has been 
destroyed on small and large length scales: at this stage spatial cooperativity is maximum, 
and one can see that the particles not only move from the borders of blocked clusters, but 
also propagate inside, breaking big clusters into smaller ones of significantly different sizes. 

We have studied the moments X{p) = (A)p. and a{p) = ^ (A^) p. — (A)p. , where A = n^/"^, 
n is the cluster size and expectation values are taken over the probability distribution P* at 
tduster- A(/9) is the typical size of the blocked clusters that slow down the dynamics, and can 
hence be seen as a (maybe rough) determination of the length S(p), defined in [10] as the 
spacing between mobile regions. A fit actually gives X{p) ~ 0.83 exp(exp(0. 12/(1 — p))). The 
fiuctuations measured by cr(/o) obey the same exponential law in p, but reach higher values 
than X{p). 

Conclusions. - We have shown that dynamical geometrical constraints in the KA model 
are able to induce spatial correlations through different characteristic length scales, that all 
diverge in the same way as p — > 1. We show in figure [S] the behaviour of ^{p), lc{p), 
and A(p). They can be rescaled by a constant, and are shown together with the fit y ~ 
exp(exp 0.16/(1 — p)). 

We have already noticed that also time scales (that grow in our allowed range far more 
than length scales, giving a far better way of discriminating between different behaviors) 
suggest a scaling law compatible with a double-exponential singularity at p = 1. We show in 
particular in figure El the time t* defined from the maximum of xit) versus p: the fit with a 
double-exponential is surely preferred over the best fit to a power divergence at pc ~ 0.881 
(note that the y scale is logarithmic, and the power fit is far from the data). 

In our study, lc{p) measures the correlation between sites of the persistence times. The 
coherence length ^(p) deduced from measurements of four-point correlation functions is the 
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Fig. 4 - For p going from 0.7 to 0.8: a) number of persistent connected clusters (lower densities on the 
left, higher densities on the right); b) probability distribution of the number of particles in persistent 
clusters at tciuster{p) for all densities: one can observe the collapse of the curves. 




Fig. 5 - The three lengths, A(p) (squares), 1.78lc{p) (diamonds), and 3.41^(p) (triangles) versus p, 
together with the best double-exponential fit. 
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Fig. 6 - t* versus p. The continuous line is for the best fit for a double-exponential singularity at 
p = 1, the dashed line is for the best power fit (with pc = 0.881). 
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typical length of cooperative events for a given density. The typical length between mobile 
particles that contribute to the relaxation of the system can be estimated via X{p), the 
average size of persistent clusters up to tduster , the time of maximum cooperativity. Moreover, 
the spatial structure contains even more complexity, since the distribution of this length obeys 
a power-law, which still remains to be understood. An analytical solution of the KA model is 
lacking, and would give a more precise understanding of the relaxation laws, the persistence 
properties, the existence of a dynamical exponent z, and the distributions of cluster sizes, 
probably to be linked with an analogy with percolation problems [6, 10, 11]. 

Finally, we want to stress that the coherence length ^(p) is a priori independent of any 
growing length-scale that is believed to exist in a system without time-translational invariance 
dynamics (aging). In the version of the KA model studied here, there is no aging due to the 
flat energy landscape. However, an aging version would be able to show how ^{p,tw) evolves 
with the age tw In particular it would be very interesting to study its behavior and that of 
t*{p,tm) in this case and compare it to experimental data in aging jammed systems [23]. 
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